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ABSTRACT 

It has been proposed that some hot molecular cores (HMCs) harbor a young 
embedded massive star, which heats an infalling envelope and accretes mass at a rate 
high enough to "choke off" an incipient HII region. This class of HMCs would mark 
the youngest phase known of massive star formation. In order to test this hypothesis, 
we model this type of object calculating the radiative transfer through a spherically 
symmetric dusty envelope infalling onto a central OB star, with accretion rates from 
M = 6 X 10~^ to 10~^ M0 yr~^. The dust thermal spectrum from infrared to radio 
wavelengths is derived and is compared with the observed fluxes of several hot cores 
which may be internally heated. We find that the data are best fitted using an envelope 
with the density distribution resulting from the collapse of a singular logatropic sphere, 
instead of that of a singular isothermal sphere. We conclude that several of these 
sources may be undergoing an intense accretion phase and find in all the cases that the 
accretion luminosity exceeds the stellar luminosity. We discuss the implications of this 
phase on the formation of massive stars. 



Subject headings: stars: formation — dust — HII regions — ISM: individual (G34.24 
+ 0.13MM, W3(H20), Orion hot core, IRAS 23385+6053) — radiative transfer 



1. INTRODUCTION 

One of the most interesting results of the past few years in the study of massive star formation 
is the discovery of Hot Molecular Cores (HMCs) which are small (^ 0.1 pc), dense {n ~ 10^ - 
10^ cm~^), hot (T ~ 100 - 300 K) and dark {Ay 10^ mag) molecular clumps. They are generally 
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found in the proximity of ultra-compact HII (UCHII) regions, and are probably directly related to 

the formation of high mass stars. The best known example is the Orion-KL hot core (e.g., Genzel 
& Stutzki 1989; Kaufman, Hollenbach & Tielens 1998 and references therein) but several other 
cores have been observed with qualitatively similar characteristics: the core near G31.41+0.31 
(Cesaroni et al. 1994b), G9.62+0.19F (Hofner et al. 1996), G34.24+0.13MM (Hunter et al. 
1998, hereafter H98), W3(H20) (Turner & Welch 1984), and IRAS 23385+6053 (Molinari et al. 
1998; hereafter M98) among others. These sources have not been detected as free-free emitters 
at centimeter wavelengths; rather, their presence is established, and their physical properties 
determined through millimeter-wave and molecular line observations (see the review of Kurtz et 
al. 1999). 

It is not yet understood what heats the gas in HMCs to a temperature higher than 100 
K. In some cases, the heating source might be a nearby star, facing the clump. This seems to 
be the case of the warm ammonia clumps associated with the source G61.48-I-0.09 (Gomez et 
al. 1995). Alternatively, Walmsley (1995) proposed that the heating source could be a recently 
formed OB-type star (or stars) inside the core undergoing an intense accretion phase. Following 
Yorke (1984), he suggested that high mass accretion rates of the infalling material could quench 
the development of an UCHII region, so that the free-free emission from the ionized material 
would be undetectable at centimeter wavelengths. Observationally, Cesaroni et al. (1994a) found 
that several HMCs seen in NH3(4,4) are coincident with groups of water masers, which are known 
signposts of massive star formation. Thus, these sources are of particular interest because they 
may represent the youngest phase yet observed in the life of a massive star, and they could help 
us understand the process of high-mass star formation. 

Kaufman et al. (1998) modeled the temperature distributions of internally and externally 
heated molecular cores and calculated the column density of hot dust and gas inside the cores. 
They considered constant density distributions and collapsing density distributions (p oc r~^/^) for 
the cores. In the latter case, they did not discuss explicitly the effect of the accretion luminosity in 
the core heating. They applied their models to the Orion hot core with mass Mcove ^ 15 Mq, and 
size Tcore ^ 0.01 pc (i.e. a column density of N{H2) > 5 x 10^^ cm~^). These authors concluded 
that an internal energy source is necessary to heat up the entire column density of the core to 
the observed temperatures, T > 200 K, since the nearest known source (Radio Source "I") is not 
luminous enough to externally illuminate and produce the observed high temperatures of the gas 
and dust. 

In this paper, we calculate the intrinsic dust thermal spectrum from infrared to radio 
wavelengths for a simple theoretical construct of the temperature profile of a HMC: a central 
massive star undergoing spherical accretion of a free falling envelope of gas and dust. The central 
heating source has two components, the stellar luminosity and the accretion luminosity. 

We compare the spectrum of our model with the observed fluxes in several sources to 
determine the physical size of the HMC, the mass accretion rate of the envelope, and the spectral 
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type of the central star. We present in §2 the parameters and the assumptions made in our model, 
in §3 the model results, in §4 a comparison with the available observations, in §5 a discussion our 
results, and in §6 the conclusions. 

2. THE MODEL 

Inspired by the scenario proposed by Walmsley (1995), we model a HMC as an envelope of 
gas and dust that is freely falling onto a recently formed massive central star. The young star is 
embedded within the dense core and interacts with it through its radiation, heating the core from 
inside. 

We assume that the system consists of a central heating source, with a total bolometric 
luminosity Lcoie, surrounded by a spherically symmetric envelope, with inner radius and outer 
radius -Rcore- The inner radius of the envelope is set by the dust destruction radius (see §2.4 ). 
The luminosity of the heating source, Lcorc = i^* + i^aco has two components, the luminosity 
of the central star (considered to be a ZAMS massive star), L*, and the accretion luminosity, 
-^acc = GM^M / R^, where G is the gravitational constant, is the stellar mass, M is the mass 
accretion rate, and i?* is the stellar radius. 

The effect of an energetic stellar wind is not taken into account. In spherical symmetry this 
is justified as long as the mass accretion rate is larger that the wind mass loss rate. Mace > Myj, 
since the wind terminal speed is of the order of the free-fall velocity. In this case the ram pressure 
of the accretion flow will prevent the wind from escaping from the stellar surface. More likely, the 
angular momentum of the infalling envelope will deviate the flow from spherical accretion and the 
material will be deposited in a circumstcllar disk around the star. In this case, the stellar winds 
could escape more easily through the poles. Our models will still be valid if the stellar winds are 
ejected in narrow bipolar cones, coexisting with the accreting envelope as in the case of low mass 
proto-stars (Adams, Lada &; Shu 1987), as long as these cones of missing material represent only 
a small fraction of the total core material and provided the axis of the bipolar outflow is not close 
to the line-of-sight. Also, the deviations of the envelope density profile from spherical symmetry 
will be important only within the centrifugal radius (Adams & Shu 1986). This region does not 
make a significant contribution to the emergent flux as long as it is small compared to the size of 
the core. Since this is the case expected for our models, the predicted spectra will be appropriate. 

To calculate the emergent dust thermal spectrum, we ignore the interaction of radiation from 
the central object with envelope matter located inside Rd- We calculate the flux received from 
the source by an observer at a distance D by integrating the equation of radiative transfer along 
rays for different impact parameters (e.g. Mihalas 1978). We take the volume emissivity of the 

dust to equal the LTE value pKi,B^{T), where p is the matter density, is the monochromatic 
absorption opacity per unit mass, Bi, is the Planck function, and T is the local dust temperature. 
In order to solve the transfer equation, a knowledge of p, and T is needed. 
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2.1. Density Distribution 

Two different types of density distributions are considered: one resulting from the collapse 
of the singular isothermal sphere (SIS; Shu 1977), and the other from the collapse of the 
singular logatropic sphere (SLS; McLaughlin Sz Pudritz 1997). The logatropic equation of state, 

P = PqIu^p/ pq), has been invoked to explain the lincwidth-sizc relation observed in molecular 
clouds (Lizano &: Shu 1989; Myers & Fuller 1992; McLaughlin & Pudritz 1996) since the sound 
speed, Cs, given by = dP/dp = Po/p, behaves like the observed velocity dispersion in these 
clouds. In the above equations, Pq has a constant value and po is an arbitrary reference density. 

The logatropic and isothermal collapses occur in the same general fashion: an expansion 
wave moves outward into a cloud at rest and the gas behind the wave falls into the central 
proto-star. The infalling matter close to the center has free-fall density and velocity profiles given 
hy p = M(327r^GM*)~V2 7.-3/2^ ^^^^ ^ ^ -(2GM*)V^r~^/^, respectively, while in the outer region 
the density profile tends towards the hydrostatic equilibrium configuration: p{r) = a^(27rG)~^r~^ 
for the SIS, where a is the isothermal sound speed, and p{r) = {Po/27rG)^^'^r^^ for the SLS. 
Nevertheless, there arc important differences in the behavior of both types of collapsing clouds. 
The mass accretion rate is constant for the SIS collapse: M* = niQa^/G, where the reduced mass 
is rriQ = 0.975. For the SLS collapse, however, the mass accretion rate is a steep function of time 
given by 

M = ^ mot^ (1) 

and the reduced mass is niQ = 0.0302. 

For a given stellar mass, M^,, and present mass accretion rate, M, the age of the system in 
the SIS is given by t^ge = M^/M, while it is longer for the SLS 

■Ml 
M 



iage = 4^ = Kse- (2) 



Furthermore, for a given central star, the mass of the infalling envelope is given by 

Menv = M,f^-lV (3) 



mo 

where m{x) is the dimensionless mass function of the dimensionless variable x. For the SIS, 
X = r/at, where r is the distance to the central star, while x = 4:(2TTGPo)~^^^rt~^ , for the SLS. For 
fage) the location of the expansion wave. Tew, occurs at x = 1. Then, one has the infalling density 
distribution within rew and the hydrostatic density distribution outside this radius. In particular, 
at the expansion wave, m(l) = 2 for the SIS and m(l) = 1 for the SLS. Thus, at a given time, 
within the expansion wave, the SLS has only 3 % of the mass in the star and the rest is contained 
in the envelope. In contrast, for the SIS, rrig = 0.975, so about half of the mass is in the star and 
half in the envelope. Also, eqns. (2) and (3) imply that, for a given and M, the isothermal 
envelope has less mass than the logatropic envelope. For simplicity, we assume that the density 
distribution is not affected by effects of the finite core size. 
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2.2. Dust Opacity 

We assume that the dust in the core is a mixture of 4 types of grains: graphite, sihcates, iron 
and water ice, with optical constants and abundances taken from D'Alessio (1996; see references 
therein), and a MRN-type size distribution (Mathis, Rumpl & Nordsieck 1977). Water ice only 
contributes to the opacity for T < 100 K. For higher temperature we assume it has sTiblimated 
(Sandford & Allamandola 1993). Since little is known about grain opacities in the sub mm range, 
we adopt a power law of the form, kx oc A^'^, with 1 < /3 < 2 for A > 200 /im. With these 
monochromatic opacities, the Rosseland mean opacity, xr, and the Planck mean opacity, np, 
are calculated in the wavelength range 0.1 < (A///m) < 10^. In this range, the monochromatic 
opacities for the low and high temperatures in the envelope (~ 20 - 2000 K), are well represented. 

2.3. Temperature Profile of the Dust Envelope 

Kenyon, Calvet & Hartmann (1993) made detailed models of the temperature structure and 
emission of dusty infalling envelopes around low mass proto-stars. They proposed a simplified 
calculation of the temperature profile which we have adopted to calculate the temperature 
distribution in our model. The basic idea of this approximate treatment considers radiative 
equilibrium for the outer optically thin dusty envelope and assumes the standard diffusion 
approximation in the inner optically thick region. The radius that divides the optically thin and 
thick regions is called the "photospheric radius", i?ph. The total luminosity is conserved at this 
radius, thus, Lcorc = 4'/ri?pjj(TTpj^, where a is Stefan-Boltzmann constant and Tph = T{Rpi^). Wc 
also require that T^{Tpi^) = 2/3, where rR(Tpii) is the Rosseland mean optical depth weighted by 
the Planck function evaluated at Tp^. Both, the optical depth and the luminosity conservation 
conditions are used to determine Rph and Tph. Unlike the SIS collapse density distribution, in 
the collapse of a SLS the photospheric radius changes for different external radii of the envelope 
because the external part of the SLS envelope has an appreciable optical depth. 

The equation for radiative equilibrium is 

/ K^B^[T{r)]du = / K^J^du. (4) 
Jo Jo 

For the optically thin outer region, the mean intensity, Juir), is approximated as 
J^(r) = W{r)B4Tph), where W{r) is the dilution factor: W{r) = i(l - [1 - (Rph/rf]^). 
Thus, one accounts for the geometric dilution of the radiation field of the envelope's photosphere 
but not the exponential attenuation of this emission which, by assumption, is less than e~^/^. 
Therefore, we consider that the grains in this region are heated only by the energy from the 
photosphere and wc neglect the heating by the diffuse field from the optically thin region. Using 
the Planck mean opacities, the thermal balance equation for the dust becomes 



(5) 
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where the Planck mean opacity, Kp{T), is evaluated at the local dust temperature, T, and Kp(rph) 
is evaluated at Tph- Once the photospheric temperature is known, the implicit eqn. (5) determines 
the temperature distribution, T{r) in the outer region. 

The innermost regions of the accreting dust envelope (inside the photospheric radius) are 
optically thick. In such regions the temperature gradient is determined by the standard diffusion 
approximation 

For a given density p(r) and Rosseland mean opacity XRiT), one solves this equation for the 
temperature distribution T{r), using a Runge-Kutta integration technique subject to the boundary 
condition T{Rp^) = Tph- The approximation used here to obtain the temperature distribution is a 
very simplified approach since a full angular and frequency-dependent radiative transfer treatment 
is beyond the scope of this paper. 



2.4. Dust Destruction Radius 

The dust is sublimated for temperatures higher than 1200 K (e.g. Adams &; Shu 1985), 
producing an inner dust free cavity. The dust destruction radius is calculated following the 
evolution of an average sized dust particle (a^ = 0.1 /xm) until it disappears, neglecting particle 
accretion onto the grains. The equation for the grain radius variation is given by 

where pg is the grain mass density, is the mean molecular weight of the gas and dust material, 
mn is the hydrogen mass and Pyap is the vapor pressure, taken from Lamy (1974). Pvap depends on 
the dust grain composition, which we have assumed to be pure silicate. With this approximation 
one obtains r(i?d) ~ 1200 K, where the dust destruction radius, R^, is of order tens of AUs (see 
Tables la, lb and 3). In a more rigorous analysis one should calculate the destruction radius of 
each one of the components and sizes of the dust grain mixture. 



2.5. Central Stars 

To determine the total luminosity, we must specify not only the stellar mass and mass 
accretion rate but also the radius and luminosity of the central star. It is unclear whether a star 
formed under very high accretion rates, will have a normal ZAMS radius. In fact, Adams & Shu 
(1985), using the results of Stabler, Shu Sz Taam (1980) for the formation of low mass stars, found 
that the stellar radius is an increasing function of accretion rate: i?* oc M°-^^. Furthermore, 
Stabler, Palla &; Salpeter (1986) studying the formation of primordial massive stars, found the 
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same type of behavior: i?* oc M°-^^. In addition, Beech & Mitalas (1994) and Bernasconi & 

Macdcr (1996) have studied the formation of massive stars with solar abundances and under 
constant accretion rates. These authors do not discuss enlarged radii, maybe because they studied 
only accretion rates < W~^Mq yr~^. It is then quite probable that OB stars, formed under 
intense accretion flows inside HMCs, will have radii larger than the ZAMS values since they have 
not had time to get rid of excess internal energy. We consider two possibilities: stars with ZAMS 
radii of Thompson (1984) (see also Vacca, Garmany, &; Shull 1996), and stars with stellar radii 
-R* = 10"*^^ cm. This last value for the radius is a factor of ~ 2-5 larger than the ZAMS radii for 
the central B0-B3 stars considered in our models and is within the range of radii found by the 
authors mentioned above. Finally, we adopt the ZAMS luminosities of Thompson (1984; see also 
Vacca, Garmany, &; shull 1996), even though, the luminosities of massive stars forming under 
intense accretion rates are also uncertain. 

3. MODEL RESULTS 

This section presents the results obtained for a set of models with different values of the 
mass accretion rate, the external radius of the envelope, and the spectral type of a ZAMS central 
star. We discuss the effects in the spectrum, when the above parameters are changed, assuming 
the collapse density distribution of the SLS and SIS. For both cases, a distance of 4.9 kpc was 
assumed for the model source, namely the distance to IRAS 23385+6053, which is one of the 
studied HMCs. 

The summary of the physical properties of the SLS models is given in Table la and its general 
trends are shown in Figure la. The top panel of Figure la shows the spectra obtained for three 
different values of the mass accretion rate (M = 1.5 x 10~^, 4.4 x 10~^ and 8.1 x 10~^ Mq yr""^), 
for a BO central star with a radius i?,^ = 4 x 10^^ cm (Thompson 1984), and an envelope with an 
external radius i2corc = 0.1 pc. As seen in the figure, the peak flux density increases and shifts to 
longer wavelengths as the accretion rate increases. This behavior can be explained as a result of 
the increase of the optical depth in the envelope, because of the overall density increase for higher 
accretion rates. As a consequence of the increase in opacity, the flux tends to be redistributed 
towards longer wavelengths, where the envelope is optically thinner, producing a shift of the peak 
of the spectrum. This also has the effect of increasing the depth of the silicate absorption feature 
at 10 fim. The total luminosity, ^core, increases with the mass accretion rate because of the 
increase in the accretion luminosity. This results in higher values of the peak Atix density. Tabic la 
shows that the photospheric radius increases and its temperature decreases with increasing mass 
accretion rates. Again this is easily understood as a result of the overall density increase in the 
envelope that moves the photospheric radius to outer (and thus cooler) regions, as the accretion 
rate increases. Table la also indicates that the dust destruction radius, i?d, increases with the 
accretion rate, since the central luminosity, i^corej increases. 

The middle panel of Figure la shows the spectra for different spectral types of the central 
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star (07, BO and B3), for an envelope with external radius of 0.1 pc, and an accretion rate of 
4.4 X 10~^ Mq yr~^. The peak flux density of the spectrum shifts to longer wavelengths for later 
spectral types. This happens because, for the models presented in this figure, i2corc < ^cw (see 
column 7 in Table la), thus, the mass of the envelope is given by eqn. (3), where m{x) oc x^'^^, 
and X oc M~^M^/^. Hence, the envelope mass is inversely proportional to the stellar mass, i.e. 
Menv oc M~°-^^ (see last column in Table la). Therefore, the opacity of the envelope increases for 
later spectral types and the emission is redistributed towards longer wavelengths. 

Finally, the bottom panel of Figure la illustrates the trends of the spectra for different 
external radii of the envelope (-Rcore = 0.15, 0.10 and 0.05 pc), for a central star with spectral type 
BO and an accretion rate of 4.4 x 10~^ Mq yr~^. As seen in the figure, the flux at long wavelengths 
increases slightly, while it decreases slightly at short wavelengths, as the external radius increases. 
This effect can be also understood as a consequence of the redistribution of the flux density to 
longer wavelengths, because of the increase of the optical depth for larger external radii. Note 
that, the photospheric radius increase with increasing external radii, as discussed in § 2.3. 

The summary of the physical properties of the SIS models is given in Table lb and its 
general trends are shown in Fig lb. The SIS models in each panel were computed with the same 
parameters that the SLS models (see columns 1 and 2 in Tables la and lb). As one can see from 
Figures la and lb, that the general behavior of the SIS models in each panel is the same the 
SLS models. However, one can also see that the spectra resulting of the SIS models have smaller 
millimeter fluxes and greater infrared fluxes than the SLS. This happens because the mass of the 
SIS envelope, Menv = M{18GM^y^/^RUre, is always smaller than that of the SLS envelope (see 
last columns in Tables la and lb). Thus, the SIS envelope is optically thinner and emits most of 
the radiation at shorter wavelengths than the SLS envelopes. 

In summary, the spectrum at all frequencies is very sensitive to the mass accretion rate; 
infrared flux densities are sensitive to the spectral type of the central star, and the value of the 
external radius has only a small effect on the spectrum for the parameters explored in Figures la 
and lb. 



4. COMPARISON OF THE MODEL WITH OBSERVATIONS 

To date, about 20 candidate HMCs have been proposed (Kurtz et al. 1999). In general, these 
objects are associated with nearby or embedded FIR sources and/or UCHII regions. Thus, it is 
very likely that the observed fluxes are contaminated by these nearby sources and overestimate 
the actual contribution of the HMC, making it difficult to fit the intrinsic spectrum. For this 
reason, we selected relatively isolated objects, and among them, those with the most complete 
data sets available in order to constrain, as much as possible, the parameters of our model. With 
these restrictions, we chose the sources: G34.24+0.13MM, W3(H20), Orion hot core and IRAS 
23385+6053 as a representative sample to test the model. 
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In order to match the sources, fits using both the SLS and SIS density distributions have 
been considered. Figure 2 shows the observed fluxes of the source G34.24+0.13MM (H98) 
and the results of tlircc model spectra for HMCs with i?corc = 0.07 pc and a central B3 
star with 10 Mq. The dotted line corresponds to the SIS collapse density distribution with 
M = 1.0 X 10~^ Mq yr~^, while the dot-dashed line corresponds to the SLS collapse density 
distribution with M = 6.5 x 10~^ Mq yr~^. For both models we used a ZAMS stellar radius, 
i?* = 2.0 X 10^^ cm (Thompson 1984). In order to match the observed fluxes at mm wavelengths, 
massive envelopes with high accretion rates are required. In particular, for a given mass accretion 
rate, the envelope is less massive for the SIS than the SLS, so one requires a higher mass accretion 
rate, producing a flux density at 20 /nm several orders of magnitude higher than the upper limit 
on the flux at this wavelength (see discussion in § 3). The logatropic collapse model, on the other 
hand, has a smaller excess emission at 20 fim which can be attenuated invoking a high external 
extinction to the HMC, Ay ~ 100 mag. This value of the external extinction is somewhat high, 
even though HMCs are found inside dense regions in molecular clouds. More generally accepted 
values are Ay ^ 30 mag (e.g. Testi et al 1998). 

The solid line in Figure 2 shows the spectrum of the SLS model discussed above, but for a star 
of radius i?* = 10"^^ cm, as discussed in § 2.5. Clearly, the excess of emission at 20 fim is reduced, 
and now flts the constraint at 20 /um, with no external extinction required, since the accretion 
luminosity decreases for this large stellar radius. Note that if this large stellar radius were used 
for the SIS density distribution, one still requires Ay > 100 mag to attenuate the excess emission 
at 20 fim. 

A similar result is found for the remaining sources, so, for these, we present only models for 
a SLS collapse distribution. In addition, the central stars are assumed to have radii i?* = 10^^ 
cm, and any external extinction is neglected. It is likely that in reality, a combination of some 
external extinction {Ay ^ 30 mag) plus a star with a radius larger than its ZAMS value, could 
be responsible for the observed spectrum. Finally, similar trends, as those shown in Figure 1, are 
found for these type of central stars with large radii. 

4.1. G34.24 + 0.13MM 

The source G34.24 + 0.13MM discovered by Hunter et al. (1998), is located 84" (1.5 pc) 
southeast of the UCHII region G34.26 + 0.15, within the G34.3 + 0.2 HII region complex, which 
is at a distance of 3.7 kpc (Kuchar & Bania 1994). H98 identified G34.24 + 0.13 MM as a massive 
proto-stellar object, which coincides with a methanol maser. This source has upper limits on the 
centimeter continuum emission of S^^ ^ 0.8 mJy at 15 GHz, and S^ ^ 0.6 mJy at 8.3 GHz (Goss 
1999, personal communication). From observations at 1.3 mm and 2.7 mm, H98 measured an 
angular diameter of 2" (0.037 pc). They also obtained near- infrared images and photometry at 
10 and 20 fim, but no infrared source was detected coincident with the millimeter position. H98 
estimated a total luminosity in the range of 1600 - 6300 Lq. These authors argue that because 
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of the high luminosity and the lack of cm compact radio continuum emission the core probably 
contains a deeply embedded proto-B star, since there is no nearby source that could provide the 
observed heating. 

We fitted the observed spectrum of this source with a model requiring a central B3 star, 
M = 6.5 X 10^*^ Mq jt^^, Rcoic = 0.07 pc, and an opacity index /3 = 2 (sec Tabic 2). Earlier 
spectral types than B3 are too luminous and are incompatible with the observed mid-infrared 
emission (see trends of the models in the middle-panel of Figure la). Accretion rates smaller 
than the value M = 6.5 x 10~^ Mq yv~^ do not reproduce the observed millimeter emission (see 
the trends of the models in the top-panel of Figure la). The top-left panel of Figure 3 shows 
the spectrum. Table 3 lists the physical properties of the source derived with this model: the 
dust destruction front, i?^; the photospheric radius, -Rph; the photospheric temperature, Tp^; the 
molecular volume density at the external radius, n(i?corc); the total luminosity, Lcore; the fraction 
of the accretion luminosity, to the total luminosity, L^cd Lcore\ the core mass within Rcore, Afcore; 
the total column density N{H2); the present age of the core given by eqn. (2); and the expansion 
wave radius. 

The required central B3 star is late enough so that it will not develop an HII region detectable 
from 3.7 kpc. Even if the star were earlier than B3, one would not expect centimeter radio 
emission since the critical mass accretion rate that would quench the development of an UCHII 
region (see Walmsley 1995, eqn. [1]) is only Merit = 3.20 x 10~^ Mq yr~"^, considerably smaller 
than that implied by our model (see Table 2). 

4.2. W3(H20) 

About 6" (0.06 pc) east of the UCHII region W3(0H) located at a distance of 2.2 kpc 
(Humphreys 1978), a luminous embedded young star was discovered near a complex of H2O 
masers by Turner &; Welch (1984; hereafter the TW object or W3(H20)). This object is manifest 
by a warm, dense concentration of molecular material (Turner & Welch 1984; Wilson, Gaume & 
Johnston 1993; Wink et al. 1994), a bipolar outflow traced by VLBI proper motions of H2O maser 
spots (Alcolea et al. 1992) emanating from a small, elongated synchrotron emission source (Reid 
et al. 1995; Wilner, Reid k Menten 1999). Additionally, Wilner, Welch & Forster (1995), and 
Wyrowski et al. (1997) detected a compact (~ 2" = 0.02 pc) mm continuum source, that was 
interpreted as thermal dust continuum emission from the W3(H20) hot core. Recent molecular 
line and continuum interferometric observations at 1.4 mm, reaching subarcsecond angular 
resolution (Wyrowski et al. 1999), indicate that the continuum emission from W3(H20) splits into 
three components. The eastern component (component A in the nomenclature of Wyrowski et al. 
1999) coincides with the source of the maser outflow and synchrotron emission, and is associated 
with a 200 K temperature peak, derived from molecular line transitions. Nevertheless, since our 
models are spherically symmetric, we model the spectrum of the compact source as observed by 
Wilner et al. (1995) and Wyrowski et al. (1997) with lower spatial resolution. 
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We fit the spectrum of this source with a model requiring a central BO star, 
M = 1.2 X 10"^ Mq yr-\ i?core = 0.05 pc, and an opacity index P = 1.6 (sec Table 2). 
Table 3 lists the physical properties of this model and the resulting spectrum is shown in the 
upper right panel of Figure 3. Because of the few data points available for this source, the model 
fit was not unique, and we selected the model solution with the spectral type of the central star 
and the index of the opacity law in accordance with those inferred by Wyrowski et al. (1997, 
1999). High angular resolution observations at mid- and far- infrared wavelengths are necessary to 
constrain the total luminosity of this source, and the spectral type of the central star. 

4.3. Orion hot core 

Because of its proximity (D = 480 pc; Genzel et al. 1981), the most studied region of 
high-mass star formation is the Orion molecular cloud (OMC-1). Line studies of the central 
portion of OMC-1 (i.e., the Kleinmann-Low nebula) showed that it was separated into a number 
of velocity components that were subsequently associated with several spatial components as 
very high-resolution mapping became available (see Wright et al. 1996 and references therein). 
Prominent at the center of this region is a dense condensation of gas and dust known as the 
"Orion hot core" which has a gas density n(H2) ^ lO'^ cm^'^ and a temperature of 150-300 K (see 
Kaufman et al. 1998 for references). Whether this source is heated either externally (Blake et al. 
1996, Wright, Plambeck & Wilner 1996) or internally (Masson & Mundy 1988) is still a matter 
of great debate. In addition. Chandler &; Wood (1997) suggested that there is an HII region 
associated with Radio Source "I" , located to the northwest of the hot core ( 2" away) . This source 
has been proposed to be the heating source of the Orion hot core. Nevertheless, as discussed in 
§1, Kaufman et al. (1998) modeled the physical conditions in the Orion hot core, and concluded 
that an internal energy source is more likely, since the Radio Source "I" is not luminous enough 
to externally heat the core up to the observed temperatures. Although high resolution molecular 
line and continuum studies have found evidence of substructure within the core (Migenes et al. 
1989; Blake et al. 1996), we apply our spherically symmetric model as a first approximation to 
this source. 

The parameters of our best-fit model are: a mass accretion rate M = 1.1 x 10"^ Mq yr~^, a 
B3 central star, an outer radius of the envelope i?corc = 0.025 pc, and an opacity index (3 = 1.6 
(see Table 2). The mass accretion rate is determined mainly by the mm and submm data. The 
spectral type of the star is constrained by the measured flux density at 30 /im (Wynn- Williams 
et al. 1984), since an spectral type earlier than B3 would have too much flux at this wavelength. 
Furthermore, the convolved image of the model at 1.3 mm has a FWHM of ~ 3", as observed 
by Blake et al. (1996; Figure 2, top panel). The total luminosity given by our model is 
Lcore = 2.5 X 10^ Lq, being most of this due to the accretion. It is worth noting that this value 
exceeds the minimum luminosity required by the models of Kaufman et al. (1998) to explain the 
observed high temperature in the Orion hot core. 
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4.4. IRAS 23385+6053 

IRAS 23385+6053 has a quite broad data set available, spanning from centimeter to near 
infrared wavelengths (see M98). No emission is detected at centimeter wavelengths above S-cr 
level of ~ 0.5 mJy. At 3.4 mm, M98 detected a compact source with a deconvolved angular size 
of 4{'5 X 3"6 (or 0.048 pc for an assumed distance of 4.9 kpc). They estimated a bolomctric 
luminosity for this source L ~ 1.6 x 10^ Lq. They also detected a compact outflow in SiO and 
HCO"'' lines that is centered on the millimeter source. 

In order to explain the lack of a detectable HII region associated with this source, M98 
consider two possibilities: a heavily obscured BO star with residual accretion, or a proto-star 
undergoing very intense accretion with M = 10~^Mq yr^^. For the latter case, using the 
relation found by Stabler et al. (1986), discussed in §2.5, they proposed a central proto-star with 
M* = 39 Mq and = 70 Rq, and assumed that the observed luminosity is dominated by the 
accretion luminosity. 

Our best-fit (model I) to the observed spectrum is shown in the lower right panel of Figure 
3 and upper left panel of Figure 4. This model I (whose parameters are given in Table 2) has a 
central B2 star, an accretion rate of 10^^ Mq yr^^, an external radius -Rcore = 0.1 pc, and an 
opacity index /?= 1.6. The physical properties of this model are given in Table 3. The lower left 
panel of Figure 4 shows the intensity profile at 3.4 mm predicted by our model and the observed 
intensity profiles along the major and minor axes of the 3.4 mm map given in Figure 1 of M98 
(dot-dashed line). As can be seen in the figure, the agreement is very good. The intensity profile 
of the model was obtained by convolving (in two dimensions) the intensity distributions at 3.4 
mm given by the model with a Gaussian beam of FWHM = 3''8 (corresponding to the size of the 
beam in the observations of M98). Note that no fit was done to the observed intensity profile at 
3.4 mm; only the observed flux densities (that are integrated values, with no information on the 
spatial distribution) were fitted with a model spectrum (upper panel of Fig. 4) . It is remarkable 
that the predicted 3.4 mm model intensity profile, obtained from the parameter set that provide 
the best-fit to the spectrum, is in such a good agreement with the observed intensity profile (lower 
panel of Fig. 4) . 

The right panels of Figure 4 shows the predictions for an alternative HMC model, similar to 
the one proposed by M98 (model II), namely a core with a central proto-star with M^, = 39 Mq, 
i?* = 70 Rq, and Lcorc = -^^acc- This model has M = 1.6 x IO^^Mq yr~^, Rcorc = 0.19 pc and 
[3 = 2. Even though the observed spectrum can be fitted by this model, the predicted intensity 
profile at 3.4 mm is too broad and its peak intensity is too low to be detectable above the rms noise 
of the map (see lower right panel in Fig. 4). This difference is due the large value of -Rcore j with 
respect to that of Model I. Therefore, this model appears to be in conflict with the observations. 
We conclude that a B2 star is most likely embedded in the IRAS 23385+6053 core, as indicated 
by model I (Table 2 and 3 and left panel of Fig 4) . 
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4.5. Temperature Profiles of HMC models 

Figure 5 shows the temperature profiles for the models presented in Figure 3 (from 
to i?core)- The different sources are labeled in each panel. These temperature profiles are 
discontinuous at the photospheric radius, i?ph, because the dilution factor, W{Rp}i) = 1/2, and, 
in our approximation, eqn. (5) then gives T < Tph at this radius. Instead, there should be a 
smooth transition between the optically thin and the optically thick regimes and hence a smooth 
temperature profile. With our simplified treatment for the temperature profile, the difference 
between the total luminosity of the spectra of model and the input luminosity, Lcore; is less than 
15 - 20 %. It can be seen that the temperature distributions in the optically thick region have 
a steeper dependence with radius than in the optically thin region, as discussed by Adams Sz 
Shu (1985) and Kenyon ct al. (1993). Finally, for all cases, at the external radius, -Rcorej the 
temperatures are T > 30K. 

5. DISCUSSION 

5.1. General Discussion 

The model studied here, of an envelope infalling onto a central massive star, reproduces 
very well the currently available observational data for the sources selected in our sample. We 

have explored two possible density distributions for the structure of the infalling envelope and 
concluded that a SLS collapse density distribution provides a better agreement with the available 
observational data, than a SIS density distribution. In order to fit the observations of these 
sources, our model requires massive, early B-type central stars (10-16 Mq), with high values of the 
mass accretion rate (^ 10~^ Mq yr~^), and core radii ^0.1 pc. The values of the mass accretion 
rate obtained exceed by a large factor the "critical" value required to quench the development of 
a detectable UCHII region, which for an early B-typc central star (the proposed central stars), is 
^ 10^^ A/q yr^^. Therefore, the model lacks a detectable HII region. Furthermore, very small 
fiux densities from dust thermal emission arc predicted at centimeters wavelengths. 

The values that we have obtained for the radius of the cores are roughly consistent with those 
inferred from the observations. However, a detailed comparison is not yet possible, since most of 
the observations give only upper limits for the source sizes. Only in the case of the 3.4 mm map 
of IRAS 23385+605 (M98) , was it possible to compare the intensity distribution predicted by the 
model with that observed, and we found a remarkably good agreement with the observed total 
flux, source size and peak intensity. Additional maps of good S /N and high angular resolution 
are required for a larger number of sources and wavelengths. In this manner, continuum intensity 
profiles may be able to discriminate between possible models. Similarly, high angular resolution 
FIR images are needed. Our model predicts a peak flux density > 10^ Jy at ~ 100 /xm. Thus, 
infrared fluxes around this wavelength, obtained with higher angular resolution than currently 
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available to avoid contamination with nearby objects, would be specially relevant in order to 
further constrain the parameters of the modeled sources. 

The intrinsic infrared colors [100-60], [60-25] and [25-12] predicted by our models (see Table 
4) are very different from those of observed UCHII regions, [100-60] =0.26, [60-25] =0.87, and 
[25-12]=0.91 (Wood & Churchwcll 1989). In particular, a color criterion [60-25] > 1.3 could be use 
to find isolated HMCs, when their fluxes are not contaminated by nearby sources. 

Furthermore, for all the sources, the accretion luminosity clearly exceed, the stellar luminosity 
(see Table 3). This result points out the relevance of the energy released in the accretion process 
at this early evolution stage, a basic ingredient that should not be ignored when modeling this 
kind of object. 



5.2. Radiation Pressure Onto Dust Grains 

The radiation pressure onto dust grains imposes an upper limit to the amount of mass that 
can accrete onto a massive star. It is therefore important to assess whether or not this effect can 
halt the collapse of the envelopes that we consider in our models. Kahn (1974) found that the 
radiation pressure of the diffuse infrared radiation field inside a dusty envelope around a massive 
proto-star decelerates the flow allowing a maximum luminosity to mass ratio L/M^ = 5436 Lq/Mq, 
although this number depends on the properties of the dust grains. Wolfire & Cassinelli (1987; 
hereafter WC87) revised the calculations of the maximum mass that an accreting proto-star can 
accumulate taking into account updated properties of the interstellar grains. Jijina & Adams 
(1996) reexamined this problem including the effect of angular momentum. In this section, we 
ensure that our models of hot cores fulfill the general constraint that guarantee that inflow will 
proceed towards the central proto-star. 

At the outer boundary of the inflow, the outward radiative acceleration must be less than the 
inward gravitational acceleration of the gas. This condition can be written as (WC87, eqn. [10]) 

pr r 



47rcGM(i?core) 

where is the flux mean radiation pressure coefficient, c is the speed of the light, M(J?core) is 
the total mass inside i?core (star plus core mass) , and -Lcore = + L^cc is the total luminosity 
available to heat and push the core material. Our models have -Lcore ^ 6 x 10^ Lq, typical core 
masses M(i?corc) ~ 130 Mq, and we approximate the flux mean opacity with the dust Planck 
mean opacity evaluated at the photospheric temperature. The photospheric temperatures we find 
are Tph ^ 100 therefore, Kp ^ 2 cm^ g~^. For this range of values 

^■^^ Ucm2g-iA6xlO4L0A 130Mo j ^ 
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and the fluid at the outer boundary can flow inward. 

On the other side, at the dust destruction front, i?d, where the direct steflar radiation is 
absorbed, the momentum rate of the flow must be greater or equal to the pressure of the steUar 
radiation fleld 

M > , (10) 

where is the velocity of the flow at the dust destruction front, R^. This condition is satisfled for 
our model (M* ~ 10 - 16 Mq, M ~ 6 x 10"^ - 10"^ Mq yr'^) as can be seen in Fig. 5 of WC87. 
Their figure also shows that our cores have less than the Eddington luminosity, LEdd, where the 
accretion is halted by radiation pressure on free electrons with an opacity given by Thompson 
scattering. 

Considering an optically thick dusty cocoon plus an optically thin envelope, and assuming 
an analytic form for the Rosscland mean opacity XR = XRo{T /Tg)'^ ■, Kahn (1974) obtained the 
maximum luminosity to mass ratio L/M^ = 5436 Lq/Mq, mentioned above. His eqn. [65] implies 
that this ratio is inversely proportional both to the cube of the dust sublimation temperature, 
Tgyj^, and to Xroj the Rosscland mean opacity for radiation at a color temperature = 22,000 K. 
Kahn used T^uh = 3675 K and XRO = 600 cm^ g^^. We find that the sublimation temperature 
for the grain mixture in our models (see §2.4) is Tg^h ~ 1200 K. Also, the value of XRO must be 
increased by a factor of ~ 8 to agree with our value for the Rosscland mean opacity at T ~ 1000 
K. These modifications, due to the currently accepted properties of the standard grain mixture 
used in this work, imply that the maximum luminosity to mass ratio derived by Kahn would 
currently be L/M* ~ 20,000 Lq/Mq. This result contrasts with the detailed modeling of the dust 
destruction process carried out by WC87. They found that very massive stars can only be formed 
with modifications of the standard MRN grain mixture for accretion onto massive stars with = 
60, 100 and 200 Mq, in which the luminosity to mass ratios are larger than L/M* = 8974 Lq/Mq. 
On the other hand, Jijina Sz Adams (1996) found that the maximum luminosity to mass ratio is 
less restrictive for a rotating collapse than for the spherical infall. In any case, all our models have 
luminosity to mass ratios Lcore/-^* < 4300 Lq/Mq, which implies that accretion can continue 
onto the central star. 

We therefore consider that radiation pressure onto dust grains will not stop the accretion flow 
onto the central proto-star for the modeled HMCs discussed in this work. For simplicity, wc have 
not considered any deviations of the flow from the dynamics of a logatropic collapsing core. 



5.3. End of the Accretion Phase 



The luminosity to mass ratios of the modeled HMCs discussed in this work are high but still 
remain below the critical ratios for radiation pressure onto dust grains needed to halt the accretion 
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(see §5.2). Nevertheless, this ratio, given by 

-^core _ GMif 

will increase with time both because the luminosity of the central star increases with time, and 
because the mass accretion rate increases as M oc (eqn. [1]) in the logatropic model. Therefore, 
radiation pressure could eventually halt the collapse of matter onto the central star. This is only 
a transient process since, as the accretion is shut off, so is the accretion luminosity, an important 
contributor to the total luminosity. Thus, one is tempted to speculate that, if a stellar wind can 
turn on after the flow is impulsively reversed, the wind could help clear out the accreting flow. An 
HII region can then be produced, although in this scenario its evolution would be governed by the 
reversal of the accretion flow, a complex problem worth studying in detail. 

If this process occurs on a short time scale compared to the accretion time scale, one can make 
a simple estimate of the final masses of the central stars inside the HMCs when the accretion flow 
is halted by radiation pressure. We assume a main sequence mass-luminosity relation, oc M^, 
and mass-radius relation, oc M^-^^ (e.g. Bowers & Deeming 1984), even though, as discussed in 
§ 2.5, it is not clear if these relations hold for massive stars formed under intense accretion rates 
which increase steeply with time, considered in this work. Figure 6 shows the luminosity to mass 
ratio as function of the dimensionless time t/tage, where tage is the present age of the HMC (solid 
line). We have assumed a 12.6 Mq central star and a mass accretion rate M = 1 x 10"^ Mq yT~^, 
at the dimensionless time t/tage = 1- These values correspond to the model for IRAS 23385+6053 
which has tagc = 5 x 10^ yr (see Tables 2 and 3). The dot-dashed line is the stellar contribution 
to the luminosity to mass ratio and the dashed line is the accretion contribution to this ratio. 
The triangles show the mass of the central star (M* oc t^, see eqn. [1]) labeled on the right 
hand vertical axis. Independently of the exact critical value for the luminosity to mass ratio for 
radiation pressure to halt the accretion flow (see discussion in §5.2), this figure shows that the 
evolution of this ratio is very fast so one would expect that in less than 10^ yr the accretion of 
mass onto the central star of the HMC should stop and an ultracompact HII region should develop 
inside these cores. 

Finally, if all OB stars arc formed within HMCs under high accretion rates, one can estimate 
the expected number of HMCs in this phase. Taking the rate of formation of O stars in the 
Galaxy N ~ 1/300 yr"^ (Wood & Churchwell 1989), and the time in the HMC phase when the 
central star is a late B star, At ~ 10^ yr, one should expect to find ^ 300 HMCs at the present 
time. On the other hand, the period of time that a core contains an early O star is very small. 
At ^2x 10* yr, so less than 60 HMCs with very luminous massive stars inside should be expected. 
The number of expected HMCs can increase by a factor of 10 if a higher star formation rate of O 
stars, ~ 4 X 10~^ yr~^ (Giisten k Mezger 1983), is used. 
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6. CONCLUSIONS 



Our main results can be summarized as follows: 



1. The observed spectra of several HMCs (G34.24 + 0.13MM, W3(H20), the Orion hot core and 
IRAS 23385+6053) can be reproduced by massive envelopes accreting onto young massive 
central stars. In order to fit the available data, we require early B-type stars with high mass 
accretion rates, M ^ 6 x 10~'^ Mq yr~^, and time-scales tage ^ 10^ yr- The main heating 
agent is the accretion luminosity (Lace > L*; see Table 3). 

2. The spectra are strongly dependent on the assumed envelope density profile. Isothermal 
collapsing envelopes have less mass in the envelope for a given mass accretion rate, with 
respect to the logatropic collapsing envelopes. Therefore, they require very high accretion 
rates (> IO'^Mq yr ^) in order to reproduce the observed millimeter emission. Nevertheless, 
high accretion rate SIS models predict strong mid-infrared emission which is, in general, not 
observed. 

3. Models of the structure of massive stars formed under intense mass accretion rates that 
increase with time, discussed in this work, are necessary to determine the characteristics 
(e.g. radius, intrinsic luminosity) of this type of object. In particular, we have used central 
stars with radii larger than the ZAMS radii in order to avoid excess fluxes at mid infrared 
wavelengths with respect to those observed. 

4. High angular resolution mid and far IR observations (which avoid contamination by nearby 
massive stars) are crucial to constrain the total luminosity of the HMCs. 

5. The critical mass accretion rate required to quench the development of an UCHII region is 
very small compared to the accretion rates required by all the HMC models studied here. 
Therefore, at this stage the HMCs have no detectable HII region. 

6. We speculate that when a critical luminosity to mass ratio is achieved, the infall will be 
halted by a powerful stellar wind and a detectable UCHII region can appear. 
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Table la. Physical Properties of a Grid of SLS Hot Core Models 



Model 


Lcore 


Rd 


Rph 




Tph 




Mcore 




(lO-^Lo) 


(AU) 


(lO^AU) 


(lO^cm-^) 


(K) 


(lO^AU) 


(Mq) 


BO star; Kcorc=0.1 pc 
















M {Mq yr-i) 
















1.5x10-3 


16 


118 


11.0 


14 


74 


21 


495 


4.4x10--* 


6.4 


47 


3.2 


2.6 


110 


47 


86 


8.1x10-5 


3.2 


18 


0.4 


0.2 


260 


144 


8 


Rcore=0.1 pc; M = 4.4 X IQ-'' Mq yr-1 
















Spectral type 
















07 


15.2 


57 


2.3 


1.1 


162 


95 


36 


BO 


6.4 


47 


3.2 


2.6 


110 


47 


86 


B3 


4.3 


45 


4.2 


4.3 


87 


30 


146 


BO star; M = 4.4 x IQ-" M© yr-l 
















Rmve (pc) 
















0.15 


6.4 


47 


3.9 


1.8 


100 


47 


205 


0.10 


6.4 


47 


3.2 


2.6 


110 


47 


86 


0.05 


6.4 


47 


2.3 


4.3 


130 


47 


19 



Table lb. Physical Properties of a Grid of SIS Hot Core Models 



Model 


Lcore 


Rd 


Rph 


n(-Rcore) 


Tph 


Tew 


Mcore 




(lO^Lo) 


(AU) 


(lO^AU) 


(lO^cm-S) 


(K) 


(lO^AU) 


(Mo) 


BO star; Rcoro=0.1 pc 
















M (Mq yr-^) 
















1.5x10-3 


16 


115 


6.5 


2.3 


98 


4.1 


141 


4.4x10--* 


6.4 


47 


2.1 


1.0 


136 


9.4 


47 


8.1x10-5 


3.2 


18 


0.3 


0.1 


284 


29 


4 


Rcore=0.1 pc; M = 4.4 x 10-"' Mq yr-1 
















Spectral type 
















07 


15.2 


56 


1.6 


0.9 


192 


19.0 


20 


BO 


6.4 


47 


2.1 


1.0 


136 


9.4 


47 


B3 


4.3 


44 


2.5 


1.0 


114 


6.0 


57 


BO star; M = 4.4 x IQ-"' Mq yr-^ 
















Rcore (pc) 
















0.15 


6.4 


47 


2.3 


0.4 


130 


9.3 


84 


0.10 


6.4 


47 


2.1 


1.0 


136 


9.4 


47 


0.05 


6.4 


47 


1.6 


4.0 


157 


9.4 


11 
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Table 2. Model Parameters Determined from Spectral Fits to Individual Sources 



Source 






M 




-Score 


/3 






(Mo) 


(Lo) 


(Mq yr- 




(pc) 




(kpc) 


G34.24+0.13MM 


10.0 


1.0x10^ 


6.5x10- 


-4 


0.07 


2.0 


3.7 


W3(H20) 


15.8 


2.5x10^ 


1.2x10" 


-3 


0.05 


1.6 


2.2 


Orion Hot Core 


10.0 


1.0x10^ 


1.1x10" 


-3 


0.02 


1.6 


0.5 


IRAS 23385+6053 


12.6 


2.8x10^ 


1.0x10" 


-3 


0.10 


1.6 


4.9 



* i?* = 10 cm is assumed. 

^ Taken from Thompson (1984), for a given M^=. 

^ The adopted distances arc taken from Kuchar &; Bania 1994, Humphreys 
1978, Genzel et al. 1981, and Mohnari et al. 1998. 



Table 3. Physical Properties of Best-Fit Models to Individual HMCs 



Source 


Rd 
(AU) 


fiph 

(103 AU) 


(K) 


rt(-Rcore) 

(10'^ cm-3) 


-f/core 

(104 Le) 


-^^acc/'^'core 
(%) 


-^core 
(Mo) 


N(H2) 
(1024 cm-2) 


^age 

(104 yr) 


(103 AU) 


G34.24+0.13MM 


35 


3.8 


72 


1.0 


1.6 


94 


131 


34 


6.1 


20 


W3(H20) 


74 


5.3 


87 


2.0 


6.7 


62 


83 


48 


5.3 


24 


Orion hot core 


52 


3.2 


88 


5.0 


2.5 


96 


27 


55 


3.7 


16 


IRAS 23385+6053 


54 


7.0 


62 


1.1 


3.1 


91 


382 


47 


5.0 


21 
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Table 4. Infrared Colors of the Best-Fit Models 

Source [100 - 60] [60 - 25] [25 - 12] 

G34.24+0.13MM 0.04 2.32 4.60 

W3(H20) -0.18 1.36 3.58 

Orion hot core -0.16 1.27 3.43 

IRAS 23385+6053 0.20 2.86 5.84 
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Fig. 1. — (a) Spectra for different sets of the SLS collapse density distribution HMC models. The 
top panel shows the trends of the spectra obtained for three different values of the mass accretion 
rate, for a BO central star and an external radius -Rcore = 0.1 pc. The middle panel shows models 
with different spectral types of the central star, for a given external radius i?corc = 0.1 pc and 
accretion rate of M = 4.4 x 10~^ Mq yr^-*^. The bottom panel shows models with different external 
radii of the envelope, for a central BO star and an accretion rate M = 4.4 x 10~^ Mq yr~^. (b) 
Same as 1(a) but for the SIS collapse density distribution HMC models. In all distance of 

4.9 kpc were adopted. 

Fig. 2.— Observed flux densities of the source G34.24+0.13MM (Hunter et al. 1998; Goss 1999, 
personal communication). Three modeled spectra with i?corc = 0-07 pc and a central B3 star 
(10 Mq) are also shown. The dotted line corresponds to the collapse density distribution of the 
SIS, with M = 1.0 X 10~^ Mq yr~^, and the dot-dashed line corresponds to the collapse density 
distribution of the SLS, with M = 6.5 x 10~^ Mq yr~\ when ZAMS radii are adopted. The solid 
line shows the spectrum of the SLS model but with a star that has a large radius, i?* = 10^^ cm. 

Fig. 3.— Observed flux densities for the sources G34.24+0.13MM (Hunter ct al. 1998; Goss 1999, 
personal communication), W3(H20) (Wyrowski et al. 1997; Wilner et al. 1995), Orion hot core 
(Wright et al. 1992; Murata et al. 1991; Masson et al. 1985; Mezger et al. 1990; The fluxes and 
errorbars at 30 /xm and 1.3 mm was estimated from the maps of Wynn- Williams et al. 1984 and 
Blake et al. 1996, respectively), and IRAS 23385+6053 (Molinari et al. 1998; IRAS-PSC2), and 
model spectra. The parameters of each model are listed in Table 2. 

Fig. 4. — Observed flux densities (Molinari et al. 1998; IRAS-PSC2) and two models for the source 
IRAS 23385+6053. The upper left panel shows the spectrum for a model I, with B2 central star, 
M = 10"^ Mq yr~-^, -Rcore = 0.1 pc, and (3 = 1.6. The lower left panel shows the intensity profile 
of the 3.4 mm emission given by the model (solid line) and the observed intensity profiles along the 
major and minor axes of the 3.4 mm map given in Fig. 1 of M98 (dot-dashed lines). The upper 
right and lower right panels correspond to model II with M = 1.6 x 10~^ M0yr~^, i?corc = 0.19 pc, 
and /? = 2, and a central proto-star with = 39 Mq and i?* = 70 Rq (for this model, Lcore = -^acc 
and is neglected) 

Fig. 5. — Temperature profiles in the range i?d < r < -Rcore; for the models shown in Figure 3. 
The different sources are labeled in each panel. The parameters of each model are given in Table 2. 
The dust destruction radius, R^, the photospheric radius, -Rph, and the photospheric temperature, 
Tph, for each source are given in Table 3. 

Fig. 6. — Luminosity to mass ratio as function of the dimensionless time t/tage (solid line). The dot- 
dashed line is the stellar contribution to this ratio and the dashed line is the accretion contribution. 
The triangles show the mass of the central star, labeled on the right hand vertical axis. 




10 100 1000 104 




10 100 1000 10* 105 




10 100 1000 10* 105 



Frequency (GHz) 




10 100 1000 104 




10 100 1000 10* 105 




10 100 1000 10* 105 



Frequency (GHz) 




10 100 1000 104 

Frequency (GHz) 




10 100 1000 10^ 10^ 10 100 1000 10^ 10^ 




